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I.  INTRODUCTION 

In  order  to  guide  ultra-intense  (>  1018  W/cm2)  laser  pulses,  plasma  channels  have  to  be 
employed.  Such  guiding  is  useful  for  applications  such  as  laser  wakefield  accelerators  [1], 
or  x-ray  lasers  [2],  The  guiding  characteristics  of  plasma  channels  are  well  known  [3-6]. 
When  focusing  the  laser  radiation  into  the  channel,  one  would  typically  like  to  choose  the 
lens  such  that  the  waist  is  located  at  the  channel  entrance,  and  such  that  the  waist  radius 
is  matched  to  the  channel  [7].  There  are  cases,  however,  where  this  cannot  be  arranged, 
and  the  question  arises  as  to  whether,  given  a  certain  focusing  configuration,  there  will  any 
guiding  at  all.  This  report  gives  a  formula  for  the  smallest  f-number  lens  that  may  be  used 
such  that  guiding  is  achieved.  It  also  extends  the  well  known  theory  of  spot  size  oscillations 
in  a  channel  to  the  case  where  the  channel  radius  is  finite. 

The  outline  of  this  report  is  as  follows.  The  finite  radius  plasma  channel  and  critical  f- 
number  are  introduced  in  section  II.  A  ray  tracing  analysis  is  employed  to  obtain  the  critical 
f-number  in  section  III.  An  eigenmode  expansion  method  is  employed  to  obtain  the  critical 
f-number  in  section  IV.  An  example  is  worked  out  in  section  V. 


II.  PLASMA  CHANNEL 


An  ideal  plasma  channel  has  a  radial  density  in  the  form 


n(r )  =  n0  (  1  + 


K, 


(i) 


where  Rch  is  the  density  doubling  radius.  Because  the  density  increases  in  r  without  bound, 
the  channel  will  guide  radiation  no  matter  what  f-number  is  used  to  couple  it  into  the 
channel.  When  the  plasma  channel  is  finite,  however,  there  is  a  critical  f-number  below 
which  the  channel  will  not  guide.  In  the  following,  we  consider  a  density  profile  of  the  form 


n{r )  = 


(2) 


n0  (1  +  r2/R2ch)  ,  r  <  Rm 
0  ,  r  >  Rm 

where  Rrn  is  the  maximum  channel  radius.  The  index  of  refraction  is  therefore  determined 
by 


v2  = 


1  -  (k2p/k2)  (1  +  r2/ R2h)  ,  r  <  Rr 
1  ,  r  >  Rr 


(3) 


where  kp  is  the  plasma  wavenumber  evaluated  at  n0  and  k  is  the  radiation  wavenumber. 
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III.  RAY  TRACING  ANALYSIS 


The  equations  for  the  evolution  of  a  ray  in  a  material  with  a  radially  varying  index  of 
refraction  are 

rtf)  rtn 

(4) 


dd  dr / 

—  =  c—  cosu 
dt  dr 


dr  ■  a 
—  =  c  sm  9 

dt 


(5) 


where  6  is  the  angle  between  the  ray  trajectory  and  the  z-axis,  and  we  assumed  r/  ~  1. 
Substituting  Eq.  (3)  for  gives 

d6 


kl 


r 


(6) 


provided  r  <  Rrn.  Taking  9  - Cl,  differentiating  once,  and  substituting  in  Eq.  (5)  gives 


where 


The  solution  to  this  equation  is 


k2  r2 

n2  =  p  c 


k2  R\h 


9  =  d-^eint  +  c.c. 


(7) 


(8) 


(9) 


It  should  be  noted  that  is  the  frequency  of  oscillation  of  a  ray,  which  is  half  the  frequency 
of  the  spot  size  oscillation.  Inserting  the  solution  for  9  into  Eq.  (5)  and  integrating  gives 


r=^eim  +  c.c. 

2iVL 


(10) 


If  r  <  Rm  for  all  t  we  can  consider  the  ray  to  be  guided.  That  is,  the  condition  for  guiding 
is  r  <  Rm,  or 


k 

9oj—Rch  ^  Rm 

kp 

Associating  angle  and  f- number  through  9  =  1/2 F  gives  the  critical  f- number  as 

1  Rch  k 

■t*  c. 


(ii) 


2  Rm  kp 


(12) 
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IV.  EIGENMODE  EXPANSION  ANALYSIS 


A  source  dependent  expansion  (SDE)  [8]  analysis  for  the  case  where  Rrn  — >  oo  has  been 
carried  out  in  Ref.  [9]  where  the  equation  for  the  spot  size  (neglecting  the  ponderomotive 
and  relativistic  terms)  was  given  as 


d2rq 


+  /  dx(l-x)  1  +  7^rX)e 


o-X 


(13) 


dz2  k2,rl  jo 

where  x  —  2 r2/r2.  To  generalize  this  equation  to  the  present  case,  we  replace  the  upper 
limit  of  the  integration  with  2 R^/r2.  Carrying  out  the  integral  gives 

d2rq  4 


dz 2  /c2r? 


+  f(ra) 


where 


and 


f(rs)  = 


k2r2  p—2Rm/rs  1 

7 TWT  - 5 - 9{rs) - 

k2R2h  r|  r.5 


9  (r  s)  —  4R^(Rc/ 1  +  -^m)  +  2Rm'rs  + 
This  can  be  put  in  the  form 


1  (  dr « 


2  Ur  i  +'/w  =  c 


where  C  is  a  constant  and 


C(rs)  = 


k2r2 


dr'sf(r'S/ 


Carrying  out  the  integral, 


V(r.)  = 


y&Aj  +  grj  -  [2(14  +  O  + 

-k2Rih 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 


Equation  (17)  is  in  the  form  of  the  total  energy  of  a  particle  in  a  potential  well.  The  first 
term  is  the  effective  kinetic  energy  and  the  second  is  the  effective  potential.  The  position  of 
the  particle,  rs,  corresponds  to  the  spot  size  of  the  beam.  The  effective  potential  for  a  leaky 
channel  ( R,m  >  Rch)  is  shown  in  Fig.  1.  The  first  local  minimum  of  V  corresponds  to  the 
matched  spot  size  tm-  If  a  beam  is  introduced  into  the  channel  with  rs  =  t'm  and  drs/dz  =  0, 
then  the  spot  size  remains  constant.  If  a  beam  is  introduced  in  some  other  condition  its 
spot  size  will  oscillate  provided  the  energy  of  the  particle,  C ,  is  below  the  potential  barrier. 
If  the  energy  of  the  particle  exceeds  the  potential  barrier  it  escapes  the  potential  well  and 
the  beam  expands  indefinitely. 
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FIG.  1:  Effective  potential  for  k/kp  =  44,  kpRch  =  12,  and  kpRm  =  20. 


The  maximum  energy  of  a  confined  particle  is  the  height  of  the  potential  barrier,  which 
is  approximately  V(Rm).  Therefore  the  condition  for  a  confined  beam  is 


+  V (rs0)  <  V(Rm) 


Z=Zq 


where  z0  is  a  reference  point  and  rs0  =  rs(z0).  Evaluating  V(Rm )  gives 


V(R,n)  =  1 1  -y  + 


k2  R2 

np  Lm. 


+ 


R 


3  R 


2  R2ch 


2  e2R2ch 


It  often  happens  that  Rrm/Rrch  S>  1  and  k^R^  1  in  which  case 


V(Rm) 

Furthermore,  if  rs0  -C  Rm,  then 


K <  ( i  _  jl\  ^ 03K RL 


k2  R-ch 


2e2 


k2R2ch 


V(r,o) 


k2r% 


2  k2  R\ 


(20) 


(21) 


(22) 


(23) 


Now  consider  the  case  where  the  beam  is  focused  at  the  entrance  of  the  plasma  channel. 
Taking  zq  to  be  the  position  of  the  channel  entrance,  drs/dz\z=ZQ  =  0,  and  the  condition  for 
a  confined  beam  becomes 

9  lkp  r*°  <  0.3^^  (24) 


k2r% 


2  k2  R2h 


k2  R. \ 


The  first  term  on  the  left  can  be  rewritten  in  terms  of  the  f-number  of  the  focusing  giving 


m  k2  \ '  r2a 


1  do 

2  RL 


(25) 
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Finally,  noting  that  the  second  term  on  the  right  is  usually  small  compared  to  the  first,  we 
obtain 


Fc  =  0.65§^  (26) 

-Lim  ftp 

which  is  the  same  as  the  ray  tracing  estimate  except  for  a  constant  factor  of  1.3. 

The  SDE  equations  can  also  be  used  to  obtain  the  matched  spot  size  tm  and  the  frequency 
of  the  spot  size  oscillation.  To  obtain  an  analytical  result,  it  is  necessary  to  take  Rrn  — >  oo 
so  that 


Voo  =  V(rs,Rn 


oo  = 


4  R2ch 


ys 


2,„2 


+  kpf 


2k‘Rlh 


(27) 


Solving  dVoo/drg  =  0  for  rs  gives  the  matched  spot  size  as 


2  R, 


rM  = 


ch 


The  wavenumber  of  the  spot  size  oscillation  is 

d2V« 


K  = 


dr ? 


4 g 


(28) 


(29) 


This  is  just  twice  the  frequency  of  the  ray  oscillation,  as  expected.  In  the  case  of  finite  Rm , 
the  equation  dV / dr s  =  0  has  to  be  solved  numerically. 


V.  EXAMPLE 


As  an  example,  consider  a  plasma  channel  with  on-axis  electron  density  of  1018  cm-3  and 
a  laser  pulse  with  wavelength  0.8  /rm.  Take  Rch  =  60  gm  and  R,rn  =  250  /rm.  Then  the 
critical  f-number  based  on  Eq.  (26)  is  6.9  and  the  period  of  the  spot  size  oscillation  based 
on  Eq.  (29)  is  1660 /kp  (about  9  mm).  To  evaluate  the  accuracy  of  these  estimates,  we  solve 
the  full  equation  of  motion 


d2r.  dV 


+  ^  =  0 


(30) 


dz 2  drs 

The  results  are  plotted  in  Fig.  2.  For  f/6  focusing  the  beam  escapes  the  channel,  while  for 
f/7  focusing  the  beam  is  guided.  Hence,  the  critical  f-number  is  between  6  and  7,  which  is 
consistent  with  the  estimated  value.  It  is  easy  to  obtain  a  more  accurate  estimate  of  the 
critical  f-number  by  numerically  solving 


V 


V(Rm) 


(31) 
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FIG.  2:  Spot  size  as  a  function  of  propagation  distance  in  the  channel  for  f/6  focusing  (solid  curve), 
f/7  focusing  (dashed  curve)  and  f/10  focusing  (dotted  curve). 

For  the  parameters  considered  here,  this  gives  Fc  =  6.94  as  opposed  to  the  analytical  estimate 
of  Fc  =  6.86.  The  frequency  of  the  spot  size  oscillation  agrees  with  the  estimated  value  in 
the  case  of  f/10  focusing,  but  is  lower  in  the  case  of  f/7  focusing.  Evidently  the  spot  size 
oscillation  slows  down  as  the  f-number  approaches  Fc. 


VI.  SUMMARY 

Radiation  propagating  in  a  finite  radius  plasma  channel  can  be  lost  if  the  focusing  used 
to  couple  the  radiation  into  the  channel  is  too  strong.  The  critical  f-number  for  guiding  is 
to  a  good  approximation 

Fc  =  0.65^^-  (32) 

rim  ftp 

where  Rch  is  the  density  doubling  radius  and  Rrn  is  the  maximum  radius  (beyond  which  the 
density  stops  increasing). 
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